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ABSTRACT 

An  1178  nm  seeded  and  1069  nm  pumped  Raman  laser  system  where  the  second  Stokes  is  amplified  in  a  1121  nm 
resonator  defined  by  high  reflector  fiber  Bragg  gratings  (FBGs)  has  the  potential  of  producing  high  output  power  of 
narrow  linewidth  1 178  nm.  However,  1121  nm  power  leakage  out  of  the  resonator  cavity  around  the  gratings  was  found 
to  impact  the  performance  of  the  laser  and  needs  to  be  dealt  with  in  order  to  obtain  high  1178  nm  output  power  levels. 
In  order  to  address  this  problem,  the  causes  of  linewidth  broadening  must  be  understood.  A  fully  nonlinear  model  has 
been  built  which  involves  propagation  of  the  spectral  wave  shape  via  the  nonlinear  Schrodinger  equation  in  addition  to 
the  Raman  processes.  It  was  found  that  increases  in  1121  nm  cavity  power,  fiber  Bragg  grating  bandwidth,  and  the 
nonlinear  index  of  refraction  n2,  as  well  as  a  decrease  in  group  velocity  dispersion  /?,  leads  to  an  increase  in  linewidth 
broadening.  It  is  concluded  that  the  magnitude  of  linewidth  broadening  seen  experimentally  can  only  be  explained  if  the 
spectral  components  outside  the  bandwidth  of  the  FBGs  are  being  amplified.  Experimentally,  1121  nm  power  leakage 
can  be  handled  by  using  a  three  wavelength  WDM  on  either  side  of  the  rare  earth  doped  amplifier.  In  addition,  usage  of 
a  fiber  having  a  high  value  for  group  velocity  dispersion  and/or  a  low  value  for  nonlinear  index  of  refraction  n2  in 
addition  to  narrower  bandwidth  fiber  Bragg  gratings  may  help  reduce  the  amount  of  linewidth  broadening. 

Keywords:  Raman  laser,  resonator  cavity,  narrow  linewidth,  nonlinear  Schrodinger  equation 


1.  INTRODUCTION 

Current  sodium  guidestar  lasers  are  either  constructed  using  slab1,2  or  fiber  laser  technology3"8.  One  way  of  achieving 
the  desired  output  wavelength  of  589  nm  for  sodium  guidestar  laser  applications  is  through  frequency  doubling  of  1178 
nm  light.  Researchers  at  the  European  Southern  Observatory  (ESO)  have  been  successful  in  constructing  a  narrow 
linewidth  sodium  guidestar  laser  based  on  frequency  doubling  of  the  output  of  a  Raman  fiber  amplifier  operating  at  1178 
nm3"9.  The  ESO  reported  achievement  of  39  W  of  narrow  linewidth  1178  nm  or  26.5  W  of  589  nm  when  frequency 
doubled.  More  recently,  44  W  of  1  MHz  linewidth  1178  nm  continuous  wave  light  was  obtained  by  Raman 
amplification  of  a  distributed  feedback  diode  laser  in  a  variably  strained  polarization-maintaining  fiber  with  a  record- 
high  optical  efficiency  of  52%  when  pumped  with  a  linearly  polarized  1120  nm  fiber  laser10.  Because  of  the  all¬ 
polarization  maintaining  configuration,  a  polarization  extinction  ratio  of  30  dB  was  obtained  and  a  20  time’s  reduction  in 
the  effective  stimulated  Brillouin  scattering  coefficient  was  achieved.  A  24.3  W  589  nm  laser  was  created  from  this  by 
using  an  external  resonant  doubling  cavity.  Finally,  22  W  of  single-frequency  1178  nm  was  obtained  from  a  counter- 
pumped  two-stage  Raman  fiber  amplifier  using  an  acoustically  tailored  fiber11. 

In  order  to  obtain  a  high  power  Raman  laser  system,  it  is  necessary  to  avoid  stressing  the  power  limited  components.  In 
the  past,  a  concept  involving  amplification  of  the  pump  wavelength  (no  signal  present)  with  a  rare  earth  doped  amplifier 
that  had  been  spliced  to  nested  Raman  cavities  for  amplification  to  higher  orders  was  reported12  and  demonstrated  by 
Nicholson.  In  this  work,  81  W  of  output  power  at  1480  nm  was  achieved  when  162  W  of  1 1 17  nm  pump,  was  injected 
into  a  120  m  long  Raman  fiber.  Because  of  the  lack  of  seeding  with  the  desired  output  wavelength,  the  Raman  process 
was  initiated  by  noise  and  the  linewidth  of  the  output  wavelength  was  broad.  As  a  follow-on  to  this,  in  2012, 
Supradeepa  demonstrated13  104  W  of  1480  nm  with  220  W  of  pump  from  a  similar  Raman  system  that  utilized  an 
1117/1175  nm  wavelength  division  multiplexer  (WDM)  to  avoid  backward  lasing  at  1175  nm  in  the  oscillator.  Finally, 
the  concept  of  seeding  with  the  desired  output  wavelength  in  conjunction  with  an  unseeded  rare  earth  oscillator  that  had 
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been  spliced  onto  a  nested  Raman  resonator  was  presented  by  Mead14  in  the  patent  literature  for  amplification  of  1480 
nm. 

The  goals  of  the  work  in  this  paper  were  to  understand  the  limiting  factors  associated  with  a  dual  1069  nm  pump  and 
1178  nm  signal  seeded  Raman  laser  system  that  involved  the  amplification  of  the  second  Stokes  in  an  1121  nm  resonator 
cavity.  Both  experimental  work  and  modeling  was  performed. 

2.  OVERVIEW  OF  RAMAN  LASER  SYSTEM 

The  system  discussed  in  this  article  is  seeded  both  with  the  zeroth  Stokes  pump  wavelength  at  1069  nm  and  the  desired 
second  Stokes  output  wavelength  at  1178  nm,  see  Figure  1.  In  this  system,  the  1069  nm  is  first  amplified  and  is  then 
Raman  converted  to  1121  nm  in  a  cavity  defined  by  high  reflector  fiber  Bragg  gratings  centered  at  1121  nm.  The  1178 
nm  signal  is  then  amplified  as  it  passes  through  the  cavity.  The  first  configuration  of  the  system  where  both  the  1069 
and  1178  nm  signals  enter  the  system  through  a  common  WDM  upstream  from  the  ytterbium  amplifier  has  been 
described  in  detail  previously13.  A  reconfiguration  of  the  original  design  to  include  the  insertion  of  more  1178  nm 
isolation  was  deemed  necessary  because  of  a  strong  1178  nm  back  reflection  that  is  about  -25  dB  relative  to  the  output 
power  level,  Figure  1.  To  enable  a  greater  degree  of  isolation  (70  dB)  with  fewer  losses,  the  1178  nm  is  injected  into  the 
laser  downstream  from  the  ytterbium  doped  amplifier.  In  addition,  a  1069/1121/1178  (3  wavelength)  WDM  is  inserted 
on  the  downstream  side  of  the  ytterbium  doped  amplifier  to  remove  any  residual  1121  nm  propagating  upstream  that 
would  favorably  compete  with  the  1069  nm  for  amplification.  In  addition,  a  similar  WDM  is  inserted  on  the  upstream 
side  of  the  amplifier  to  remove  any  1121  nm  which  traversed  the  amplifier  to  avoid  damage  to  the  expensive  1069  nm 
source.  In  this  system,  because  the  fiber  Bragg  gratings  (FBGs)  are  high  reflectors,  the  1069  nm  will  be  Raman 
converted  to  1121  nm  which  will  grow  in  power  in  the  1121  nm  in  the  resonator.  The  linewidth  of  the  1121  nm  in  the 
cavity  will  broaden  to  at  a  minimum,  the  maximum  allowed  by  the  fiber  Bragg  gratings.  The  1178  nm  seed  upon 
entering  the  1121  nm  resonator  cavity  will  then  be  amplified  via  Raman  conversion  of  the  first  Stokes  at  1 121  nm  to  the 
second  Stokes  at  1178  nm.  Narrow  linewidth  output  is  achievable  because  the  system  was  seeded  with  the  desired 
narrow  linewidth  output  wavelength.  In  addition,  the  system  is  comprised  of  entirely  polarization  maintaining  10/125 
fiber. 


30  dB  1069  nm 


Figure  1.  Schematic  diagram  of  a  single  cavity  Raman  system  that  is  pumped  with  1069  nm  and  seeded  with  the  second  Stokes 
at  1178  nm. 


3.  EXPERIMENTAL  RESULTS 

Because  of  rapid  growth  of  the  1178  nm  output  due  to  initiation  of  the  Raman  process  with  a  seed  versus  noise  in 
addition  to  the  expected  very  high  power  levels  of  the  first  Stokes  at  1121  nm  in  the  Raman  resonator,  it  was  originally 
thought  that  a  Raman  fiber  (either  silica  or  germanosilicate)  having  a  length  on  the  order  of  5-20  m,  depending  on  the 
level  of  the  1178  nm  seed,  was  possible.  This  would  have  been  the  case  if  the  1 121nm  power  was  building  in  the  cavity 
versus  leaking  spectrally  around  the  high  reflector  fiber  Bragg  gratings.  Using  the  original  configuration  described  in 
reference  15,  several  sets  of  experimental  results  were  obtained  for  three  different  cavity  lengths  (11  m,  61  m,  and  111 
m).  The  fiber  utilized  in  the  cavity  was  Nufern  FUD-4112  which  is  a  10/125  polarization  maintaining  germanosilicate 
fiber.  For  even  very  short  cavities,  an  extraordinary  amount  of  linewidth  broadening  of  the  1121  nm  in  the  resonator 
cavity,  Figure  2,  was  observed. 
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The  net  effect  of  the  linewidth  broadening  was  a  large  amount  of  1121  nm  power  leakage  around  the  gratings  in  the 
spectral  domain.  This  resulted  in  less  buildup  of  1 121  nm  power  in  the  cavity  than  expected,  i.e.,  it  was  leaking  out  both 
ends  of  the  resonator  cavity,  a  reduction  in  the  conversion  of  1069  nm  to  1121  nm  as  it  traversed  the  resonator  cavity, 
and  a  lower  output  power  of  1 178  nm  than  expected.  A  longer  resonator  cavity  was  found  to  increase  the  conversion  of 
1069  to  1121  nm  and  was  found  to  have  slightly  less  linewidth  broadening  due  to  the  lower  power  levels  of  1 121  nm  in 


Figure  2  -  (a)  Linewidth  broadening  in  1 1 1  m  cavity  as  a  function  of  the  1069  nm  pump  power  above  the  threshold,  and  (b) 
Linewidth  broadening  as  a  function  of  the  cavity  length. 

the  cavity.  Although  a  longer  cavity  enabled  greater  Raman  amplification  of  the  1178,  the  ultimate  downside  would  be  a 
lower  threshold  for  SBS  and  ultimately,  lower  output  power  levels  of  1178  mn.  Figure  3.  In  the  configuration  shown  in 
Figure  1,  nearly  .45  W  of  1178  nm  was  obtained  as  an  output  prior  to  terminating  the  experiment  prematurely.  The  big 


0  5000  10000  15000 


P1069  nm  (entering  cavity)  [mW] 


(c) 


P1178nm  (output)  [mW] 


(d) 


Figure  3  -  (a)  Leakage  power  level  of  1 121  nm  from  a  grating,  (b)  Output  power  level  of  1 178  nm  exiting  the  cavity,  and  (c) 
Unused  1069  nm  ALL  as  a  function  of  the  power  level  of  1069  nm  entering  the  resonator  cavity,  (d)  Power  out  of  the  back  tap 
on  the  lower  power  side  of  the  amplifier  as  a  function  of  the  1178  nm  output  power. 

concern  was  the  backward  traversing  power  which  was  monitored  out  of  the  .01%  tap  coupler  on  the  lower  power  side  of 
the  rare  earth  doped  amplifier.  It  was  found  that  98.6%  of  the  backward  traveling  power  was  1178  nm  which  could  be 
detrimental  to  the  1178  nm  source  which  only  had  about  35  dB  of  isolation  protecting  it  in  this  configuration.  Fear  of 
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damage  to  the  1178  nm  source  is  what  drove  a  premature  stop  to  the  experiments  at  less  than  maximum  output  power 
levels  of  1 178  nm.  In  these  experiments,  some  linewidth  broadening  of  the  1178  nm  was  seen  with  the  111m  cavity,  see 
Figure  4.  The  amount  is  hard  to  discern  because  of  the  limited  resolution  of  the  optical  spectrum  analyzer. 


Figure  4  -  Spectrum  of  the  1 178  nm  output  at  various  power  levels.  At  the  highest  output  power  level,  about  1.66%  of  the 
1178  nm  output  was  initiated  from  noise. 

To  protect  the  1178  nm  source,  the  system  was  reconfigured  as  shown  in  Figure  1  with  70  dB  of  isolation  protecting  the 
1178  nm  source  from  the  back  reflection.  In  addition,  to  offset  the  additional  losses  arising  from  the  increased  isolation, 
the  1178  nm  source  was  placed  downstream  from  the  ytterbium  doped  amplifier  as  far  as  possible.  With  161  m  of 
germanosilicate  fiber  comprising  the  cavity,  an  1178  nm  output  of  1.1  W  was  obtained  which  was  close  to  the  threshold 
for  Stimulated  Brillouin  Scattering  when  only  one  temperature  zone  is  used.  Obviously,  to  move  forward,  it  is  necessary 
to  understand  the  driving  factors  behind  the  huge  amount  of  1121  nm  linewidth  broadening  in  the  resonator  cavity  as 
well  as  how  to  mitigate  it. 


4.  DESCRIPTION  OF  NONLINEAR  MODEL 

Previously  a  nonlinear  model  was  developed  that  propagated  the  spectral  wave  shape  in  frequency  space.  This  model 
was  found  to  be  applicable  for  shorter  resonator  cavities  and  for  the  case  where  the  amount  of  linewidth  broadening  was 
somewhat  limited.  This  was  the  situation  since  a  certain  minimum  number  of  spectral  points  were  needed  to  adequately 


Figure  5  -  Percentage  of  1 121  nm  power  leaking  past  a  3  nm  bandwidth  fiber  Bragg  grating  having  a  reflectivity  of  99.91%  for 
cavity  lengths  of  1 1 ,  6 1 ,  and  111  m. 

represent  the  computational  domain  and  the  number  of  computations  scaled  as  the  number  of  spectral  points  to  the  fourth 
power.  The  model  was  found  to  become  computationally  intractable  for  cases  involving  longer  resonator  cavities  where 
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there  is  a  lot  of  linewidth  broadening.  Even  with  its  shortcomings,  this  model  led  to  some  interesting  insights  in  that 
longer  cavities  don’t  necessarily  have  increased  leakage  of  1121  nm  power.  Experimentally,  the  dependence  on  the 
percentage  loss  of  1 121  nm  power  seemed  to  depend  on  intracavity  power  levels  rather  than  the  cavity  length.  Figure  5. 

To  enable  investigation  of  a  system  having  a  lot  of  linewidth  broadening,  it  is  necessary  to  utilize  the  nonlinear 
Schrodinger  equation  (1)  to  describe  the  propagation  of  a  spectral  envelope  down  a  fiber.  Incorporated  within  the 
nonlinear  Schrodinger  equation  are  the  following  nonlinear  effects:  dispersion,  four  wave  mixing,  self-phase  modulation, 
cross-phase  modulation,  etc.  In  the  case  of  this  model,  the  number  of  spectral  points  needed  is  only  those  necessary  to 
adequately  describe  a  continuous  envelope  in  discrete  space.  In  addition,  none  of  the  calculations  scale  as  the  number  of 
spectral  points  to  the  fourth  power. 

In  the  present  model,  the  spectral  envelope  is  considered  for  the  1121  nm  (1st  Stokes)  primarily  because  of  issues 
associated  with  the  huge  amount  of  linewidth  broadening  in  the  resonator  cavity.  The  spectral  envelope  is  also 
considered  for  the  1178  nm  (2nd  Stokes)  because  of  the  importance  of  linewidth  broadening  when  it  comes  to  certain 
applications.  The  1069  and  1242  nm  are  each  handled  as  a  single  wavelength  since  linewidth  broadening  is  not  expected 
to  be  significant  nor  is  it  of  interest.  First,  the  1121  nm  spectral  distribution  is  defined  by  a  set  of  discrete  spectral 
regions  having  power  levels  representative  of  the  level  of  noise  present  in  the  system  and  the  spectral  distribution  of  the 
1178  is  assumed  to  be  a  Lorentzian  representative  of  the  seed  with  noise  overlaid.  For  each  distribution,  the 
computational  domain  is  defined  to  be  broad  enough  so  that  effects  arising  from  the  model  running  up  against  the 
boundaries  of  the  calculation  are  minimized.  Next,  the  nonlinear  parameter  y  and  the  group  velocity  dispersion  f}2  are 
found  for  the  set  of  discrete  points  representing  the  spectral  envelopes  for  1121  and  1178  nm.  (See  references  16-19  for 
details  on  the  calculation  of  yand  J32.)  Once  the  initial  spectral  envelopes  have  been  defined,  the  spectral  envelopes  for 
the  1121  and  1178  nm  are  propagated  one  step  along  the  resonator  cavity  in  the  forward  direction  using  the  nonlinear 
Schrodinger  equation  and  new  spectral  envelopes  (1121  and  1178  nm)  are  obtained.  Next,  Raman  is  considered  for  the 
four  Stokes  orders  involved  (1069,  1121,  1178,  and  1242  nm)  for  the  same  step  and  the  power  in  the  spectral 
components  of  both  the  1121  and  1178  nm  envelopes  is  updated  based  on  the  power  and  Raman  gain  coefficients 
associated  with  each  individual  component.  The  process  just  described  is  then  repeated  until  the  calculation  reaches  the 
end  of  the  resonator  cavity.  At  this  point,  the  1121  nm  spectral  envelope  propagating  in  the  forward  direction  is 
reflected  off  of  the  fiber  Bragg  grating  and  the  backward  traveling  power  in  the  four  Stokes  orders  is  then  propagated 
backward  along  the  resonator  cavity  using  the  process  just  described.  Once  the  fiber  Bragg  grating  is  reached  on  the 
opposite  end  of  the  resonator  cavity,  the  1121  nm  spectral  envelope  propagating  in  the  backward  direction  is  once  again 
reflected  and  the  process  just  described  is  repeated  until  convergence  is  obtained.  The  initial  version  of  the  code 
assumes  that  the  forward  and  backward  propagating  waves  do  not  interact  with  each  other  spectrally  through  the 
nonlinear  Schrodinger  equation. 

The  nonlinear  Schrodinger  equation  is  defined  as  follows17' 20: 

°Ega  =  (B  +  R)E(z,t)  (1) 


where  D  is  a  differential  operator  that  accounts  for  dispersion  and  losses  within  a  linear  medium,  IV  is  a  nonlinear 
operator  that  governs  the  effect  of  fiber  nonlinearities  on  wave  propagation,  and  E  (z,  t)  is  the  complex  field  envelope  at 

step  z  and  time  t.  D  and  N  are  given  by: 


re-i  Pn 


gm 

at™ 


(2) 


and 


N  —  iy  \E(z,  t)  | 2  + 


o)0E(z,t )  dt 


(\E(z,t)\2E(z,  t) 


(3) 


where  a  is  the  absorption  coefficient,  /?,„  are  the  dispersion  parameters  of  various  orders,  and  y  is  a  nonlinear  parameter 
associated  with  the  fiber  involved.  The  angular  frequency  at  the  center  of  the  envelope  is  given  by  <u0. 
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For  a  step  along  the  resonator  cavity  in  z,  both  the  linear  and  nonlinear  parts  of  the  nonlinear  Schrodinger  equation, 
separately,  have  analytical  solutions  but  the  nonlinear  Schrodinger  equation  as  a  whole  does  not  have  a  general 
analytical  solution.  If  the  step  size  h  taken  along  z  is  small  enough,  then  the  two  parts  of  the  nonlinear  Schrodinger 
equation  can  be  treated  separately  with  only  a  “small”  numerical  error  that  is  second  order  in  the  step  size  h.  That  is, 

E{z  +  h,  t)  «  exp(h.D )  ■  exp(hiV)  ■  P(z,  t).  (4) 

For  the  nonlinear  part  of  the  nonlinear  Schrodinger  equation:  =  iy\E\2E  —  NE,  one  obtains  the  following  analytic 

solution  for  a  small  step  h  in  the  z  direction: 

En  (z  +  h,t)  =  exp  (iy  |  E  (z,  t)  | 2  h)  E  (z,  t )  =  exp  ( WV)  E(z,  t).  (5) 

This  is  considered  to  be  an  intermediate  result  that  is  not  representative  of  the  time  domain  at  the  end  of  a  step  of  length 
h.  In  order  to  obtain  a  final  result,  dispersion  must  be  considered.  The  dispersion  step  has  an  analytical  solution  in  the 
frequency  domain  so  it  is  necessary  to  first  Fast  Fourier  transform  En  (z  +  h,  t)  into  the  frequency  domain.  Also,  the 
dispersion  operator  has  the  following  form  when  represented  in  the  frequency  domain: 

Dm=HD}  =  ~-Zm=2i^Pmmm.  (6) 

Since  the  system  is  continuous  wave,  only  the  absorption  a  and  the  group  velocity  dispersion  /?2  are  considered  in 
D(ia)).  Third  order  and  higher  dispersion  parameters  are  ignored.  The  solution  to  the  linear  part  of  the  nonlinear 
Schrodinger  equation  in  the  frequency  domain:  qp'  =  D(iui)E  at  the  end  of  a  step  of  length  h  therefore  has  the  following 
form: 


E(z  +  h,  u>)  =  exp(hD(ico))T(exp(hN)E(z,t)).  (7) 

Upon  Fourier  transforming  this  result  into  the  time  domain,  a  representation  of  the  time  domain  at  the  end  of  the  step  of 
length  h  is  obtained: 


E{z  +  h,t)—T  1  p(z  +  h,  (8) 

Equations  (7)  and  (8)  are  therefore  considered  to  be  final  results  representative  of  the  frequency  and  time  domains  at  the 
end  of  a  step  of  length  h. 

After  the  1121  and  1178  nm  spectral  envelopes  have  been  propagated  a  step  of  size  h  along  the  resonator  cavity  using 
the  nonlinear  Schrodinger  equation,  the  Raman  process  is  then  applied  to  the  same  step  for  the  four  Stokes  orders 
involved  and  the  resultant  increments  in  power  are  spread  over  the  1121  and  1178  nm  spectra  at  the  end  of  the  step.  For 
the  Raman  calculation,  the  total  power  in  the  1121  and  1178  nm  spectra  is  rolled  up  into  the  center  wavelength  of  each 
spectral  distribution.  An  effective  Raman  gain  coefficient  representative  of  each  spectrum  is  also  calculated  for  the  1121 
and  1178  nm.  Once  this  has  been  accomplished,  the  equations  shown  below  for  each  Stokes  order  involved  are  used  to 
propagate  the  power  a  step  h  along  the  Raman  resonator17: 


dz 
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In  these  equations,  the  upper  sign  in  each  refers  to  the  forward  propagating  wave  and  the  lower  sign  refers  to  the 
backward  propagating  wave.  Also,  yp  refers  to  the  Raman  gain  coefficient  for  a  pump  wavelength  of p  nm.  Note  that 

yp  =  9r  1  1000  nrn  scales  inversely  with  the  pump  wavelength  and  the  gR  in  this  equation  refers  to  the  Raman  gain 

coefficient  for  a  pump  wavelength  of  1  micron.  In  addition,  P/  refers  to  the  forward/backward  propagating  power  in  the 
LP01  mode  at  s  nm,  Yp  "  refers  to  the  back-scattered  Raman  gain  coefficient  at  a  pump  wavelength  of  p  nm  which  is 
taken  to  be  equivalent  to  yp,  yB  is  the  SBS  coupling  factor,  as  is  the  absorption  coefficient  at  .v  nm  and  Ay  refers  to  the 
effective  area  for  modal  overlap  of  the  pump  and  signal  wavelengths  at  p  and  s  nm. 

Once  the  Raman  calculation  has  been  performed  and  the  increase  or  decrease  in  power  is  known  for  each  Stokes  order, 
for  the  1121  and  1178  nm,  the  increment  in  power  is  then  spread  over  the  discrete  point’s  representative  of  the  spectral 
envelope  prior  to  the  next  step  with  the  nonlinear  Schrodinger  equation.  The  Raman  gain  coefficient  at  the  center 


Grating  Reflection  Spectrum 


Figure  6  -  Reflectivity  of  Teraxion  fiber  Bragg  grating  used  in  the  experiment.  The  three  analytic  regions  on  the  long  and 
short  wavelength  sides  are  labeled. 

wavelength  of  each  discrete  spectral  region  along  with  the  power  in  each  region  determines  how  the  increment  in  power 
is  spread  across  each  spectrum. 
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Finally,  once  the  end  of  the  resonator  cavity  has  been  reached,  the  power  in  1121  nm  will  be  reflected  according  to  the 
reflectivity  profile  of  the  grating  used  in  the  experiment.  Figure  6.  The  reflectivity  at  the  other  Stokes  orders  is  assumed 
to  be  negligible  in  the  current  version  of  the  model.  The  reflectivity  profile  of  the  fiber  Bragg  grating  was  modeled  by 
fitting  both  the  long  and  short  wavelength  sides  with  three  functions,  two  second  order  polynomials  for  regions  1  and  3 
and  a  linear  function  for  region  2.  In  the  model,  the  red  points  can  be  specified  to  enable  investigation  of  the  effect  of 
various  FBG  edge  reflectivity  profiles  on  linewidth  broadening.  When  a  spectral  envelope  reflects  off  of  a  fiber  Bragg 
grating,  it  is  assumed  that  the  envelope  reflects  according  to  the  reflectivity  characteristics  of  the  FBG.  After  the  power 
is  reflected  from  a  grating,  the  power  in  the  Raman  components  is  then  reassigned  to  the  backward  propagating 
components  in  order  to  enable  propagation  of  the  four  Stokes  orders  in  the  opposite  direction.  The  process  just 
described  is  repeated  as  the  light  propagates  backward  in  the  Raman  resonator.  Upon  reflection  off  of  the  other  grating, 
the  entire  process  is  repeated  until  convergence  is  obtained. 

4.1  Results  of  the  model 

In  order  to  gain  an  understanding  of  what  impacts  linewidth  broadening,  the  effect  of  the  following  on  the  amount  of 
linewidth  broadening  associated  with  a  cavity  of  length  150  m  was  investigated  by  the  model:  1.  the  group  velocity 
dispersion  (32,  2.  the  nonlinear  index  of  refraction  n2,  3.  the  1121  nm  power  in  the  cavity,  and  4.  the  high  reflectivity 


Figure  7  -  (a)  Effect  of  the  group  velocity  dispersion  on  linewidth  broadening,  (b)  Effect  of  the  nonlinear  of  index  of 
refraction  n2  on  linewidth  broadening,  (c)  Effect  of  1 121  nm  cavity  power  on  linewidth  broadening,  (d)  Effect  of  the  grating 
bandwidth  on  linewidth  broadening  ALL  as  determined  by  the  percentage  of  1121  power  beyond  the  bandwidth  of  the 
gratings.  Note  that  the  following  parameters  were  associated  with  the  calculation  for  each  case:  (a)  100  W  1121  nm  cavity 
power,  3  nm  bandwidth  gratings,  n2  =  3.73  x  10‘20  m2/W  (the  value  for  a  10%  germanium  doped  fiber  core),  h  =  .25  m 
calculation  step  (note  that  no  wavelength  dependence  of  (32  was  considered  in  this  calculation  except  for  the  case  that  was 
specific  to  the  fiber  involved  -  square),  (b)  100  W  1 121  nm  cavity  power,  3  nm  bandwidth  gratings,  (32  =  0,  h  =  .25  m,  (c)  3  nm 
bandwidth  gratings,  (32  =  0,  n2  =  3.73  x  10‘20  m2/W,  h  =  .25  m,  and  (d)  100  W  1 121  nm  cavity  power,  (32  =  0,  n2  =  3.73  x  10"20 
m2/W,  h  =  .25  m.  The  cavity  length  was  150  m  for  all  cases. 

bandwidth  of  the  FBG,  see  Figure  7.  Since  the  spectrum  is  effectively  reset  at  each  reflection  to  that  contained  within 
the  bandwidth  of  the  FBG,  one  pass  through  the  cavity  is  representative  of  the  amount  of  linewidth  broadening  in  the 
resonator  cavity  if  the  forward  and  backward  propagating  waves  are  assumed  not  to  interact  spectrally  through  the 
nonlinear  Schrodinger  equation.  As  is  widely  known,  more  linewidth  broadening  occurs  when  the  value  of  the  group 
velocity  dispersion  is  low  versus  high.  This  is  the  case  because  the  phase  matching  condition  associated  with  four  wave 
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mixing  is  broken  when  the  value  associated  with  the  group  velocity  dispersion  high.  In  order  to  get  a  feel  for  the 
maximum  amount  of  linewidth  broadening  that  is  achievable,  the  group  velocity  dispersion  was  taken  to  be  zero  when 
investigating  the  effect  of  n2,  the  1121  nm  cavity  power,  and  the  bandwidth  of  the  grating  on  the  linewidth  broadening. 
As  can  be  seen  in  Figure  7,  the  following  tends  to  result  in  increased  linewidth  broadening:  1.  an  increase  in  n2,  i.e.,  more 
germanium  in  the  fiber  core,  2.  a  greater  1121  nm  cavity  power,  3.  a  longer  cavity  length,  4.  a  greater  bandwidth  of  the 
fiber  Bragg  grating,  and  5.  a  lower  value  of  the  group  velocity  dispersion.  Of  particular  note  is  that  when  a  realistic, 
wavelength  dependent  value  was  used  for  the  group  velocity  dispersion,  linewidth  broadening  was  significantly  reduced. 

Illustration  of  the  amount  of  linewidth  broadening  as  a  function  of  the  distance  traversed  along  a  resonator  cavity  along 
with  the  spectral  shape  of  the  power  after  traversing  a  150  m  fiber  is  shown  in  Figure  8  as  a  function  of  the  intracavity 
power  level  of  1121  nm.  As  expected,  the  amount  of  linewidth  broadening  as  a  function  of  distance  traversed  increased 
faster  than  linear  for  all  1 121  nm  intracavity  power  levels.  In  addition,  the  amount  of  spectral  broadening  also  increased 
as  the  1121  nm  intracavity  power  level  increased.  Similar  trends  are  seen  for  changes  in  n2,  the  bandwidth  of  the  FBGs 


Wavelength  [nm] 


Figure  8  -  (a)  Amount  of  1 121  nm  power  leakage  outside  the  high  reflectivity  spectral  bandwidth  of  the  FBG  as  a  function  of 
the  distance  traversed  along  the  cavity,  (b)  Spectral  composition  of  the  1121  nm  power  in  the  cavity  after  traversing  150  m 
along  the  cavity  as  a  function  of  the  1121  nm  power  in  the  cavity,  and  (c)  Spectral  composition  of  the  1121  nm  power  which 
leaks  past  a  FBG  having  a  reflectivity  of  99.91%  and  a  3  nm  bandwidth.  Note  that  (32  =  0,  n2  =  3.73  x  10"20  m2AV,  h  =  .25  m, 
and  the  bandwidth  of  the  grating  was  3  nm. 

and  the  group  velocity  dispersion.  The  bottom  line  is  that  for  the  parameters  associated  with  the  experiment,  i.e.,  cavity 
lengths  of  11  to  111m,  10%  germanium  doping  of  the  core  which  results  in  an  n2  =  3.73  x  10"20  m2/W,  wavelength 
dependent  (32  having  a  value  of  1.811  x  1026  ps2/km  at  1121  nm,  3  nm  bandwidth  FBGs,  and  intracavity  1121  nm  power 
levels  of  17-30  W  (depending  on  the  cavity  length  being  studied),  even  with  p2  taken  to  be  zero  which  is  not  physically 
real,  the  amount  of  linewidth  broadening  seen  experimentally  (up  to  16-25%  1121  power  leakage  past  each  grating 
depending  on  the  cavity  length)  can’t  be  duplicated.  The  fact  that  p2  has  a  wavelength  dependence  and  is  nonzero,  tends 


Figure  9  -  Spectrum  of  power  emerging  from  a  111  m  cavity  between  1000  and  1300  nm. 
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to  limit  the  amount  of  linewidth  broadening  on  a  single  trip  up  the  fiber  quite  dramatically.  In  addition,  experimentally, 
for  one  pass  through  the  cavity,  very  little  linewidth  broadening  was  seen  for  the  1069  nm  whereas  for  the  1121  nm 
which  is  effectively  reset  spectrally  after  each  reflection  from  a  FBG  and  therefore  effectively  travels  the  same  distance, 
this  was  not  the  case.  Figure  9.  Because  of  the  dramatic  amount  of  linewidth  broadening  seen,  something  else  is 
obviously  occurring  to  support  and  enable  amplification  of  1121  nm  spectral  components  outside  the  high  reflectivity 
bandwidth  of  the  FBG.  Some  possibilities  are:  1.  Spectral  interaction  of  1121  nm  power  propagating  in  one  direction 
with  power  propagating  in  the  other  direction  in  the  resonator  cavity,  i.e.,  the  power  reflecting  off  of  a  grating  which  has 
been  spectrally  narrowed  would  see  power  propagating  in  the  opposite  direction  which  has  been  broadened  due  to 
propagation  along  the  cavity  since  the  system  is  continuous  wave,  2.  Residual  reflectivity  of  power  out  beyond  the  high 
reflectivity  bandwidth  of  the  grating,  or  3.  A  scattering  process  at  the  fiber  Bragg  grating. 

The  effect  on  linewidth  broadening  due  to  an  increased  baseline  reflectivity  out  beyond  the  high  reflectivity  bandwidth 
of  the  fiber  Bragg  gratings  was  investigated.  Note  that  data  from  the  manufacturer,  Figure  6,  seems  to  indicate  a 
reflectivity  in  these  regions  of  -25  to  -30  dB.  The  fall  off  of  the  FBG  reflectivity  to  baseline  was  assumed  to  occur  over  a 
0.2  nm  range  on  either  side  of  the  high  reflectivity  part  of  the  grating  for  modeling  purposes.  The  following  case  was 
considered:  11m  cavity,  3  nm  FBGs,  p2  =  0,  n2  =  3.73  x  10"2l),  h  =  .5  m,  1069  nm  input  power  of  40W,  1178  nm  input 


beta  2  =  0 


Figure  10  -  Percentage  of  cavity  power  leaking  past  a  fiber  Bragg  grating  in  the  forward  direction  versus  the  percentage 
reflectivity  of  the  fiber  Bragg  grating  beyond  the  high  reflectivity  bandwidth.  The  case  considered  was:  11m  cavity,  40  W 
1069  nm,  1  W  1178  nm,  3  nm  bandwidth  FBGs,  beta_2  =  0,  n2  =  3.73  x  10'20,  h  =  .5  m,  186  nm  of  calculational  bandwidth, 
and  fall  off  of  grating  reflectivity  over  .2  nm  range  on  either  side  of  grating. 

power  of  1  W,  and  186  nm  of  calculational  bandwidth.  Multiple  passes  of  the  light  through  the  resonator  cavity  was 
simulated  until  convergence  was  achieved.  As  can  be  seen  in  Figure  10,  the  percentage  of  1121  cavity  power  residing 
outside  the  high  reflectivity  bandwidth  of  a  FBG  increased  quite  dramatically  once  the  baseline  reflectivity  of  the  FBG 
reached  greater  than  80  percent.  (Note  that  the  usage  of  p2  =  0  is  the  high  linewidth  broadening  case.)  Because  the 
baseline  reflectivity  of  the  FBGs  used  in  the  experiment  was  in  the  neighborhood  of  -20  to  -30  dB  according  to  the 
specifications  and  measurements  of  the  supplier,  the  residual  reflectivity  of  the  1121  nm  spectral  components  outside  the 
high  reflectivity  bandwidth  of  a  FBG  is  not  the  source  of  the  tremendous  amount  of  linewidth  broadening  observed 
experimentally,  i.e.,  calculations  using  residual  reflectivities  of  -20  to  -30  dB  do  not  support  this. 

5.  CONCLUSION 

The  leakage  of  1 121  nm  power  from  the  resonator  cavity  is  a  problem  that  needs  to  be  dealt  with  in  order  to  obtain  high 
output  power  levels  of  1178  nm.  Experimentally,  1121  nm  power  leakage  can  be  handled  by  using  a  three  wavelength 
WDM  on  either  side  of  the  rare  earth  doped  amplifier.  In  addition,  usage  of  a  fiber  having  a  high  value  for  the  group 
velocity  dispersion  and/or  a  low  value  for  the  nonlinear  index  of  refraction  n2  would  probably  lead  to  less  linewidth 
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broadening.  In  addition,  the  usage  of  narrower  bandwidth  FBGs  might  help  as  well.  Finally,  additional  development  of 
the  model  is  necessary  in  order  to  gain  a  full  understanding  of  the  source  of  the  linewidth  broadening  in  the  1121  nm 
resonator  cavity.  In  particular,  the  investigation  of  the  effect  of  a  spectral  interaction  between  the  forward  and  backward 
propagating  waves  in  the  resonator  cavity  is  necessary  in  conjunction  with  a  more  realistic  representation  of  the  fiber 
Bragg  gratings.  To  conclude,  1121  nm  power  leakage  from  the  resonator  cavity  will  probably  drive  the  usage  of  a 
longer  Raman  cavity  than  previously  planned  since  lower  1121  nm  cavity  power  levels  leads  to  reduced  conversion  of 
the  various  Stokes  orders.  The  impact  of  this  will  be  the  need  for  mitigation  of  Stimulated  Brillouin  Scattering  in  order 
to  raise  the  output  power  level  of  the  1178  nm. 
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